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We report a numerical study of equilibrium phase-diagrams and interfacial properties of bulk and 
confined colloid-polymer mixtures using grand canonical Monte Carlo simulations. Colloidal parti- 
cles are treated as hard spheres, while the polymer chains are described as soft repulsive spheres. 
The polymer-polymer, colloid-polymer, and wall-polymer interactions are described by density- 
dependent potentials derived by Bolhuis and Louis [Macromolecules, 35 (2002), p. I860]. We com- 
pared our results with those of the Asakura-Oosawa-Vrij model, that treats the polymers as ideal 
particles. We find that the number of polymers needed to drive the demixing transition is larger for 
the interacting polymers, and that the gas-liquid interfacial tension is smaller. When the system 
is confined between two parallel hard plates, we find capillary condensation. Compared with the 
AOV model, we find that the excluded volume interactions between the polymers suppress capillary 
condensation. In order to induce capillary condensation, smaller undersaturations and smaller plate 
separations are needed in comparison with ideal polymers. 



I. INTRODUCTION 

Mixtures of colloids and polymers [TJ |5J [3] are simple 
model systems that have been studied extensively in the 
past years. Provided the size and the number of polymers 
are sufficiently high, such mixtures can phase-separate 
into a colloidal gas phase that is poor in colloids and rich 
in polymers, and a colloidal liquid phase that is rich in 
colloids and poor in polymers. While it is well estab- 
lished that the interactions between sterically stabilized 
colloidal particles are well described by the hard sphere 
potential, ^ the interactions between the polymers in the 
mixture can be in general very complicated. Neverthe- 
less, if we consider the case of flexible polymer chains in 
a "good solvent" conditions, we can assume that the ex- 
cluded volume interactions between the chains are small 
and that the chains cannot penetrate the colloidal par- 
ticles. The mechanism behind the demixing transition 
is then easily explained. In Fig. [ija), we illustrate the 
mixture of spherical colloids and polymer chains. Around 
each colloid there is a depletion region prohibited to the 
polymers due to the hard-core interactions. If two col- 
loids approach each other, so that two depletion zones 
overlap there is an increase in free volume for the poly- 
mer chains, i.e. an increase in entropy. The increase in 
entropy can be described by an effective attractive inter- 
action between the colloidal particles. A similar depletion 
mechanism occurs between a hard wall and the colloidal 
particles. If the polymers do not adsorb at the wall they 
are excluded from a region close to the wall. The overlap 
between the depletion zone at the wall and the depletion 
zone of one colloid induces an increase in free volume 
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and, hence, in entropy. 

One particular simple model for colloid-polymer mix- 
tures is the Asakura-Oosawa-Vrij (AOV) model [H [5J |7] 
that describes the polymer chains as spherical parti- 
cles with a radius equal to the radius of gyration of 
the polymer. Furthermore, polymer spheres can freely 
overlap, while they are excluded from a centre of mass 
distance from the colloidal particles. The AOV model 
has been studied extensively in the past years, and it 
was shown that it describes qualitatively the bulk [8,9, 
EH ED H HI HI E3 EES] and interfacial phase behav- 
ior Q3 OH HH mi OS) [22] of mixtures of conoids and 
polymers. A similar level of agreement was found for the 
interfacial tension of the gas- liquid [16j [23 , 24J and wall- 
fiuid interfaces, [2H [251 125] and for the phase behavior of 
confined systems. [13 HH1 12E 1101 13H 1321 1331 IS3 

The quantitative discrepancies between experimen- 
tal results and the AOV model results can arise from 
a number of reasons, like non-ideal solvent condi- 
tions, [35 colloid-induced polymer compression, |36j ef- 
fect of charges on the colloidal surface, [37J IMj or poly- 
mer excluded volume interactions. In this article, we 
will concentrate on the latter aspect. The simplest in- 
clusion of polymer interactions was done by introduc- 
ing a step function interaction between the polymers, 
i.e. an energy penalty for the overlaps of two poly- 
mers. The step potential was used to study the bulk 
phase diagram and interfacial tension [25j [39l HQ] as well 
as the stability of the fioating liquid phase in sediment- 
ing colloid-polymer mixtures [H] with a geometry-based 
density functional theory (DFT). This approach gives re- 
sults that are in better agreement with experiments when 
compared against the AOV model, but the height of the 
step potential must be introduced as an additional free 
parameter. Furthermore, we expect the polymer-colloid 
interaction to be modified as well when considering ex- 
cluded volume interactions between polymers. Other the- 
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oretical approaches have been developed to study the ef- 
fect of polymer interactions in colloid-polymer mixtures. 
Aarts et al. [32 extended the free volume theory [9] to 
include polymer interactions, and studied the gas-liquid 
interfacial tension with the square gradient approxima- 
tion approach. [42] 133] HU 05] A one-component per- 
turbative DFT that includes excluded volume interac- 
tions that uses the approach of Ref. [3S] was developed 
by Moncho-Jorda et al. [37] to study confined systems 
and the gas-liquid tension. 

Another approach is to describe the polymers as 
soft spheres, [3H1 SSI ISH] with effective interactions de- 
rived from inversion of the centre of mass (CM) cor- 
relation functions in lattice Self-Avoiding- Walk (SAW) 
polymer simulations. This approach generates soft 
density-dependent potentials for the polymer-polymer 
and colloid-polymer interactions, that give accurate sim- 
ulation results for the bulk phase behavior, [51] that are 
in quantitative agreement with experimental results. In 
a similar approach, proposed by Jusufi et al. |52j . the 
effective potentials are derived from off-lattice molecular 
dynamics simulations of SAW polymer chains. These po- 
tentials were used in Monte Carlo simulations to study 
the bulk phase behavior [53] [53] and the gas-liquid inter- 
facial tension [55] of colloid-polymer mixtures, leading to 
results in quantitative agreement with experiments. 



II. MODEL 

The colloids are treated as hard spheres and the cor- 
responding pair potential reads 



/3v cc (Rij) 



oo if Rij < a c 
otherwise, 



(1) 



where Rij = \R\ — Rj 



is the distance between two col- 
loidal particles, with Ri the position of the centre of 
mass of colloid i, (3 = l/fc^T, ks Boltzmann's constant, 
and T the temperature. For the coarse-grained effective 
polymer-polymer, colloid-polymer, and wall-polymer po- 
tentials we use the expressions of Ref. [35] obtained from 
microscopic simulations of SAW polymer chains consist- 
ing of 500 segments on a lattice and a radius of gyration 
i? 9 =16.83 lattice units at zero concentration. We intro- 
duce the overlap concentration p* defined by the equa- 
tion 4/3irp*Rg = 1, and the polymer reservoir packing 
fraction rjL = p^/p*, with p r p the density in the reser- 
voir of pure polymers in osmotic equilibrium with the 
two-component system of interest. 

The effective density-dependent polymer-polymer in- 
teractions, HE] read 



k=l 



In this work, we study the effect of excluded volume 
interactions on the interfacial properties and phase be- 
havior of confined colloid-polymer mixtures with Monte 
Carlo simulations. We simulate a binary mixture with 
the density-dependent potentials derived in Ref. [48], and 
compare the simulation results with those for the AOV 
model. 



where R g is the radius of gyration, = \fi — fj\ is the 
distance between two polymers, with fl the position of 
the centre of mass of polymer i. The density-dependent 
parameters are linear in the density = a1 + a\rfp, and 
6fe = b° + b\rf p . All coefficients, except 6|, are given 
in Table [i] The coefficients b\ are fixed by imposing the 
equality of the mean field equation of state 



f3P/ Pp = l + p p pv(0;p p )/2 



(3) 




FIG. 1: (a) Illustration of a mixture of colloids and polymer 
chains in contact with a hard wall. Depletion zones (dashed 
lines) and overlap zones (light grey) are also depicted, (b) Il- 
lustration of a model colloid-polymer mixture in the polymer- 
as-soft-sphere approach. 



for the fitted potentials and the SAW simulations, where 
the function 



Pv(0; p p ) = 4tt 



r 2 f3v pp 



(r,p p )dr 



(4) 



is the k = component of the Fourier transform of the 
polymer-polymer pair potential. In practice, the condi- 
tion is satisfied by imposing the equality between 



em*;) 



' 3/2 E 

i=l 



r\3 



(5) 



derived from Eq. ((2) and 

(3v(0;r] r p ) = 4tt(1.2902+ 

0.28132 rf p + 0.13676 {rf p ) 2 - 0.040892 (t^) 3 ) (6) 

derived using the potentials obtained from the inversion 
of the radial distribution function determined from SAW 
simulations. 
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TABLE I: Coefficients for the density-dependent parameters 
of the polymer-polymer interaction potential v pp defined in 
Eq. 



k=l k=2 k=3 
a° k 1.47409 -0.23210 0.63897 
aj -0.07689 0.03132 0.24193 
b° k 0.98137 0.42123 
b\ -0.05681 -0.02628 



TABLE III: Coefficients for the density-dependent parame- 
ters of the wall-polymer interaction potential v wp defined in 
Eq. j8j. 





k=0 


k=l 


k=2 


k=3 


ft 


62.7242 


-6.4093 


2.5081 


-0.6904 


fk 


56.4595 


-3.8880 


5.1562 


-1.5519 


fi 


-29.9283 


2.0442 


-2.1336 


0.5973 



The concentration-dependent colloid-polymer poten- 
tial reads [5T] 



k=l 

cxp[-((l^ - fj\ - e k W P ))l{R s dk.W P )))\ (7) 



where \R\ — fj\ is the distance between colloid i and 
polymer j. The density-dependent parameters are lin- 



ear in density, i.e 
e k = eg 
size ratio q 



and 



c k = 4 + c l%> dk = dk + a k%' 
„. The coefficients are given in Table II for 
1.05. 



TABLE II: Coefficients for the density-dependent parame- 
ters of the colloid-polymer interaction potential v cp defined 
in Eq. |7| for size ratio q=1.05. 





k=l 


k=2 


cl 
cl 


5.5610 
-0.8042 


1.8477 
1.4759 


dl 

dl 


0.7751 
-0.1151 


1.2720 
0.1052 


el 
el 


0.4082 
0.1410 


0.0 
0.0 



The interaction between colloidal particles and the 
hard wall is hard-sphere-like, that is the colloidal parti- 
cles cannot penetrate the walls. The interaction between 
polymers and the hard wall [48^ reads 



exv[fx{rf p )z/R g + 



f 2 ( V ;)z 2 /Rl + f 3 ( V ;)z 3 /R%(8) 



where z is the distance between the wall and the centre 
of mass of the polymer. The parameters have a quadratic 
density-dependence f k [rQ = /£ + flrf p + fi(Vp) 2 , with 
fc=0, 1, 2, 3. The coefficients are given in Table [TTT] In 
Fig. [2] we show the effective polymer-polymer, colloid- 
polymer, and wall-polymer interactions for rf v — 1.02995 
as an example. 

A final note on the potentials described in this section 
is in order. Due to a small error in the calculation of the 
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FIG. 2: Interaction potentials for polymer density rf v = 
1.02995. (a) Pair potential between polymers fiv pp (full line), 
and between polymers and colloids (3v cp (dashed line), (b) 
Wall-polymer potential /3u wp as a function of the distance 
z/a c of the polymer centre of mass from the wall. 



radius of gyration in Ref. |1H] , all the equations and co- 
efficients presented here were parameterized assuming a 
radius of gyration i? ff =16.495 lattice units, instead of the 
correct value of i? g =16.83. The simulations we carried 
out for this work, as well as the simulations presented 
on Ref. [51] , were done imposing a size ratio of q = 1.03 
(corresponding to the radius of gyration i? 9 =16.495) in 
the calculations. All the results have been interpreted 
using the correct value for the size ratio q = 1.05 (cor- 



4 



responding to the radius of gyration _R 9 =16.83) . This 
correction only change the value of the polymer packing 
fractions r\ v and rf , while the colloid packing fraction i} c 
is unaffected. 



III. METHOD 



eta r p . On the other hand, the acceptance probabilities 
for the insertion/deletion of polymers were always larger 
then 40%. The low insertion/deletion probabilities of 
colloidal particles have only an effect on the efficiency 
of the algorithm. As shown by the smooth probability 
distributions in Fig. [3] our simulations were long enough 
to get good data. 



We carried out Monte Carlo simulations in the grand 
canonical ensemble, i.e. with fixed volume, temperature, 
and chemical potentials \x c and /i p of colloids and poly- 
mers, respectively. For each value of \i v we determined 
the potentials given in the previous section, at the poly- 
mer reservoir packing fraction r) p (p p ) calculated by in- 
version of the equation of state 



R 3 g = log(p^) 



0.04658- 



+ 35.48(p^i?^) 2 - l5.7l(p r p Rp d . (9) 

This equation was derived by integrating the Gibbs- 
Duhcm equation with the pressure given by Eq. ^ and 
Eq. (§. 

To study phase coexistence, we sample the probabil- 
ity P(N c )\ Zc „t of observing iV c colloids in a volume V 
at fixed colloid fugacity z c and fixed polymer reservoir 
packing fraction rf p , using the successive umbrella sam- 
pling [56J. We use the histogram reweighting technique 
to obtain the probability distribution for any z' c once 
P{Nc)\z c ,rr is known for a given z c : 




FIG. 3: Logarithm of the probability P(N C ) (not normalized) 
as a function of the colloid packing fraction rj c for a simu- 
lation box with dimensions 12x12x16 a 3 at varying polymer 
reservoir packing fraction 7^=1.14, 1.20, and 1.23, from top 
to bottom. All state points are at coexistence. 



lnP(iV c )| 2 



In P(N C 



AUn 



(10) 



At phase coexistence, the distribution function P(N C ) 
becomes bimodal with two separate peaks of equal area 
for the colloidal liquid and gas phases. We determine 
which z' c satisfies the equal area rule 



(N c ) roo 

P{N c )\ z ,. % dN c = P{N c )\ z , K dN c 

ii J{N C ) 



with the average number of colloids 



(11) 



The liquid-gas interfacial tension j\ g is obtained from 
P(N c )\z' c ,-ri r at coexistence [57] 



7ig = 



2L 2 



In 



^(■^cmax) ~r~ ^(-^c.max) 



-ln(P(AT c>min )) 



(13) 

) and P{N l c max ) are the maxima of the 
gas and liquid peaks, respectively, and P(N Cim - ln ) is the 
minimum between the two peaks. 



where P{N« 



N c P(N c )\ z , v rdN c 



using the histogram reweighting equation (10 1. The 



(12) 



simulations are carried out in a rectangular box V — 
L x L x H, and the sampling of the probability ratio 
P(N) / P(N + 1) is done, in each window, until the differ- 
ence between two successive samplings of the probability 
ratio is smaller than 5 x 10 -4 . An example of the sampled 
probability distributions is given in Fig. [3j 

We used single particle insertion/deletion of colloids 
and polymers. The typical acceptance probabilities for 
the insertion/deletion of colloids were between 4% to 1% 
(from low to high colloid density) for state points close 
to the critical points and from 0.1% to 0.01% at high 



IV. RESULTS 



In Sec. |TT] we explained the straightforward, but non- 
trivial procedure for generating the interaction poten- 
tials. It is therefore important to check the internal con- 
sistency of our calculations. Fig. [4] shows the predictions 
of the equation of state ^ plotted against simulations re- 
sults of a grand canonical simulation of pure polymers in- 
teracting with the potential ^ . In the range of chemical 
potentials that are relevant for the gas-liquid separation 
the simulation results are consistent with the equation of 
state. 
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FIG. 4: The chemical potential P^-pot as a function of the 
polymer packing fraction r\ v in a system of pure polymers. 
The equation of state (J9j) (solid line) is compared with the 
results of grand canonical Monte Carlo simulations (circles) 
of pure polymers interacting with the interaction potential 
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Bulk phase behavior and gas-liquid interfacial 
tension 



In Fig. [5] we present the bulk phase diagram obtained 
from grand canonical Monte Carlo simulations using suc- 
cessive umbrella sampling and histogram reweighting. In 
particular, Fig. |5ja) shows the phase diagram in the poly- 
mer packing fraction r] p , colloid packing fraction r\ c repre- 
sentation. These results are consistent with the findings 
of Ref. [ST]. The free volume theory [S5] extended to 
include excluded volume polymer interactions, |42) over- 
estimates the simulation results by almost a factor of 
two. This result may be explained by the renormalisa- 
tion group theory expression used in Ref. [42 to evaluate 
the polymer interactions, that underestimates the corre- 
lation length of the polymers. Since the polymers are ef- 
fectively smaller, a higher number of polymers is needed 
to drive the demixing transition. Also shown are the 
experimental results of de Hoog and Lekkerkerker [12"] . 
The experimental polymer concentration is much larger 
than our simulation results. This discrepancy can be ex- 
plained by considering the depletion force measurements 
of Wijting et al. [SJ5] in the same colloid-polymer system 
as was used in Ref. [12]. They found that the depletion 
forces are much smaller than expected, probably due to 
adsorption of the polymers on the colloidal surface. We 
stress that the potentials used in this work compare well 
with the experiments of Ramakrishnan et al. [60] at a 
size ratio q=0.67 [ST]. 

Fig. |5jb) shows the phase diagram in the polymer 
reservoir packing fraction rf p , colloids packing fraction r\ c 
representation. The discrepancy between our results and 
those of Ref. [ST] are due to a slightly different equation 
of state used for the conversion of the chemical poten- 
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FIG. 5: Bulk phase diagram of colloid-polymer mixtures 
with a size ratio g=1.05 for volume V=516 at (circles) and 
V=2304 a\ (squares). Also shown are the binodals of the 
AOV model with size ratio q = 1.0 (dashed lines), and the 
results of Ref. |51) (solid line) . (a) Polymer packing fraction 
rip, colloid packing fraction r\ c representation. Shown are also 
the results of the free volume theory with polymer interac- 
tions |42| (dotted-dashed line) and the experimental results 
of de Hoog and Lekkerkerker [12] (diamonds), (b) Polymer 
reservoir packing fraction rf v , colloid packing fraction r) c rep- 
resentation. 



tial /ip to the polymer packing fraction in the reservoir 
r]p. In this work, we inverted Eq. (l9|, while in Ref. 
[ST] the original SAW equation of state was used. The 
binodal has a critical point at lower 77^, and the density 
difference between the gas and liquid phases increases 
for increasing rf v . This phase diagram is equivalent to 
the temperature-density phase diagram of a simple fluid, 
with the polymer reservoir packing fraction playing the 
role of inverse temperature. 

In Fig. [6] we present the simulation results of the di- 
mcnsionless interfacial tension /3j g ia^ for the gas-liquid 
interface, as a function of the difference in packing frac- 
tions between the gas and the liquid phase. The inter- 
facial tensions decreases in the case of excluded volume 
interactions with respect to the AOV model. The com- 
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parison between our results and the experiments of Aarts 
et al. [21] is quantitatively better than the results of the 
AOV model, although de Hoog and Lekkerkerker [T2] 
show that it is difficult to obtain accurate interfacial ten- 
sion measurements. In addition, we compare our results 
with the predictions of the extended free volume the- 
ory plus a square gradient approximation to evaluate the 
tension, and the density functional theory (DFT) 
of Moncho-Jorda et al. [17]. The DFT uses effective one- 
component pair potentials between the colloids that in- 
cludes the excluded volume interaction according to the 
approach of Louis et al. [3B]. The predictions of the two 
theories are very close to each other and to our simulation 
results. 



tion is in general larger for larger wall separations, and 
consequently the polymer adsorption is smaller. There- 
fore, for larger wall separations we expect a larger col- 
loid density and a smaller polymer density. The combi- 
nation of these two effects renders the interpretation of 
the phase diagram in the r\ p and rj c representation fairly 
complicated. Fig. [7]Jb) displays the phase diagram of 
colloid-polymer mixtures in the polymer reservoir pack- 
ing fraction 77^, colloid packing fraction r\ c representation. 
The critical points of the confined systems shift towards 
higher rf v for decreasing wall separation. The interpre- 
tation of the binodals in the (f7p,f? c ) representation is 
more straightforward because only r\ c suffers of the prob- 
lems described above. We can say that the huge shift 




FIG. 6: Dimensionless interfacial tension Pja^, between the 
gas and liquid phase as a function of the difference in pack- 
ing fractions A 7/ = r\i — r\ g of the coexisting liquid (r\i) and 
gas (rig) phase. Results for the interacting polymers with size 
ratio q = 1.05 (circles) are compared with the results for the 
AOV model with size ratio q — 1.0 (square). Triangles denote 
experimental results of de Hoog and Lekkerkerker [12! (trian- 
gles up) and Aarts et al. [3T] (triangles down) . The thick con- 
tinuous line indicates the DFT predictions of Moncho-Jorda 
et al. [17], while the dashed lines are the predictions of the 
square gradient approximation theory of Aarts et al. [45] . 
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Fig. [7J shows the phase diagram of colloid-polymer 
mixtures confined between two hard walls with separa- 
tion distance H/a c = 00, 16, 8, 4, 2. In particular, 
Fig. [7ja) shows the phase diagram of colloid-polymer 
mixtures in the polymer packing fraction rj p , colloid pack- 
ing fraction r\ c representation. The binodals in Figj7^a), 
hardly change under confinement. We stress that the 
comparison between the absolute densities inside capil- 
laries of different sizes is complicated due to two factors. 
First, the density for confined systems is a ill-defined 
quantity because it depends on the definition of the vol- 
ume. Our choice of volume depends on the wall separa- 
tion. Therefore, the comparison between densities for dif- 
ferent wall separations in not entirely consistent. Second, 
the adsorption of colloids and polymers is different for dif- 
ferent wall separations. We find that the colloid adsorp- 



FIG. 7: Phase diagram of colloid-polymer mixtures confined 
between two hard walls with distance H/o c — oo(bulk), 16, 8, 
4, 2. Solid lines are a guide to the eye. (a) Polymer packing 
fraction r\ v , colloid packing fraction r) c representation. Inset: 
Blow-up of the high r\ c region of the binodal. (b) Polymer 
reservoir packing fraction rf v , colloid packing fraction r) c rep- 
resentation. 

in gas densities at coexistence for the confined systems 
with respect to the bulk binodal is due to the formation 
of liquid layers at the walls (See Fig. [8ja) and (c)). The 
adsorption of colloids at the walls is due to the depletion 
attraction, which was also observed in the AOV model. 
However, the thickness of the liquid layer was in that 
case much smaller than in the present case and the shift 
in gas density less pronounced. [27J EI] 

Fig. [9] displays the phase diagram in the polymer 
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FIG. 8: Typical configurations from computer simulations of 
the coexisting gas (a),(c) and liquid phase (b),(d) for the con- 
fined system with separation distance H/a c =16 and chemical 
potentials fi c = 10.38 and fi p = 5.6 (?7p = 1.086). Colloids are 
dark grey and polymers are light grey. For clarity (c) and (d) 
display the same configurations of (a) and (b), respectively, 
without the polymers. 




10 



Kelvin Equation 
GO Bulk 
□ H/a c =16 

O H/a=8 
A H/a=4 
* H/a =2 



11 



12 



5.75- 



5.5- 



l 1 , 



13 




A H/a =4 



chemical potential \i v , colloid chemical potential /i c rep- 
resentation. These variables do not depend on the def- 
inition of the volume and arc therefore independent of 
the wall separation distance. This is the ideal represen- 
tation to compare the binodals for different wall sepa- 
ration distances. The binodals collapse to a single line 
because of the thermodynamic equilibrium conditions of 
the gas-liquid coexistence. Regions above the binodal are 
gas-like, while regions below the binodal are liquid- like. 
We find a shift of the binodals towards higher polymer 
chemical potential and smaller colloid chemical potential 
indicating the occurrence of capillary condensation. Our 



estimates for the critical points are reported in Tab. IV 



TABLE IV: Critical values of the polymer reservoir packing 
fraction rf v , colloid packing fraction rj c , and chemical poten- 
tials /j,p and fj, c for wall separation distances H/a c = oo, 16, 
8 , 4, and 2. 



H/a c 




(/^p)cr (/-£c)cr 


oo 


1.06(5) 0.192(5) 


4.99(5) 9.23(5) 


16 


1.08(5) 0.243(5) 


5.10(5) 9.45(5) 


8 


1.09(5) 0.235(5) 


5.16(5) 9.50(5) 


4 


1.17(5) 0.219(5) 


5.72(5) 10.25(5) 


2 


1.27(5) 0.199(5) 


6.43(5) 10.64(5) 



Also shown in Fig. [9] are the predictions of the Kelvin 
equation [31] 
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FIG. 9: Phase diagram of colloid-polymer mixtures confined 
between two hard walls with distance H/a c = oo, 16, 8, 4, 2, 
in the polymer chemical potential fi p , colloid chemical poten- 
tial fi c representation. In (b) we show a blow-up of the phase 
diagram. For clarity the results of H/ a c = 2 are not shown. 
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The predictions of the Kelvin equation are in good agree- 
ment with the simulation results for H/a c =16, and 8, 
but underestimates the shift for H/a c =4, and 2. 

In order to compare the results with those for the AOV 
model we scale the binodals by the b ulk c ritical point 
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(M P )J? r ulk and (p c )c? lk , reported in Tab. _ 
and (b) show that the shift of the binodals and critical 
points is smaller for the model with interacting polymers 
than for the AOV model studied in Refs. [571 EE] for all 
state points considered. As shown in Fig. [5] the difference 
p\ — p 9 c between the colloid packing fraction of the liquid 
and the gas at coexistence is smaller for the interacting 
polymers than for the AOV model for low 77L but larger 
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for high 77^. On the other hand, Fig. 
show that the shift in chemical potential is always smaller 
for the interacting polymers than for the AOV model. 
Therefore, we deduce from Eq. (14 1 that the difference 
7wi — 7wg (liquid-wall and gas-wall interfacial tensions 
at coexistence) is smaller for the interacting polymers 
model. 

Moncho-Jorda et al. H7J have shown that at constant 
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rfp the excluded volume interactions increase the differ- 
ence in wall tensions. This is not in contrast with our 
findings, because for capillary condensation we need to 
consider the difference in wall tension at bulk coexistence, 
that is at different rf v for the AOV model and the inter- 
acting polymer model. 




FIG. 10: Phase diagram of colloid-polymer mixtures confined 
between two hard walls in the polymer chemical potential 
A 1 p/(A i p)cr llk ) colloid chemical potential At c /(/i c )J?" lk represen- 
tation, (a) H/a c — oo, 4, and 2 for the model discussed in 
this paper, (b) H/a c — oo, 5, and 2 for the AOV model. 
Results are taken from Ref. [31| . 



V. CONCLUSIONS 

We have investigated bulk and confined colloid- 
polymer mixtures, using Monte Carlo simulations. Col- 
loids are treated as hard spheres, while polymers were 
described as soft repulsive spheres. Colloid-polymer, 
polymer-polymer, and wall-polymer density-dependent 
interactions were described by the coarse-grained poten- 
tials derived in Ref. [48]. We find a bulk phase behavior 
consistent with the findings of Bolhuis et al. |51j . Our re- 
sults for the bulk phase behavior are also similar to those 
for the AOV model with interacting polymers [42], but 
here the binodal line lies at higher polymer packing frac- 
tions, i.e., the number of polymers needed for the demix- 
ing transition is larger. These results arc in agreement 



with the findings of other authors [HI 051 [5TJ [55] . The 
comparison of our phase diagram with experiments |12j 
is found to be poor for the size ratio q = 1.05 considered 
here. This is surprising, since the same interaction po- 
tentials provided good agreement with experiments |60] 
at a smaller size ratio [51j . In fact, this discrepancy can 
be explained by considering the results of Wijting et al. 
|59j on depletion potential measurements on the same 
colloid-polymer mixtures that was used in the phase be- 
havior experiments. These measurements concluded that 
the depletion attraction was smaller than expected, prob- 
ably due to polymer adsorption on the surface of the col- 
loids. 

On the other hand, better agreement is found for the 
gas-liquid intcrfacial tension when compared to the ex- 
periments of Aarts et al. [21] for the same system. Our 
results show that the gas-liquid interfacial tension is 
smaller for the interacting polymers than for the AOV 
model. This is in agreement with the works of oth- 
ers on colloid-polymer mixtures with interacting poly- 
mers [101 05J 071 [55J . Both the square gradient approx- 
imation and the DFT provide a good description of the 
simulation results. 

In addition, we studied the phase behavior of the mix- 
ture confined between two parallel hard walls with sep- 
aration distance H/a c =16, 8, 4, and 2. We find that 
the hard walls induce capillary condensation, and that 
the theoretical predictions of the Kelvin equation are 
in reasonable agreement with the simulation results for 
H/a c =16 and 8, but underestimate the binodal shifts 
for H/a c =4 and 2. Compared to the AOV model the 
excluded volume interactions suppress the capillary con- 
densation. This implies that the effect can only be ob- 
served at statepoints close to bulk coexistence and that 
smaller plate separations are needed to induce capillary 
condensation for fixed supersaturation compared with 
non-interacting polymers. In other words the Kelvin 
length is smaller for interacting polymers than for the 
AOV model. In addition, we observed the formation of 
rather thick wetting layers at the largest wall separation 
we studied. This effect seems to be enhanced by the pres- 
ence of the excluded volume interactions. At large wall 
separations the wetting layers provide an effective con- 
finement H — 2t, with t the thickness of the wetting layer, 
as was recognized by Derjaguin |61j . which increases ef- 
fectively the Kelvin length [30 for the colloid-polymer 
mixtures with excluded volume interactions. The in- 
fluence of excluded volume interactions on the wetting 
properties of colloid-polymer mixtures is currently under 
investigation. 
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